How harmonic is dipole resonance of metal clusters? 
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Abstract 



We discuss the degree of anharmonicity of dipole plasmon resonances in 
metal clusters. We employ the time-dependent variational principle and show 
that the relative shift of the second phonon scales as iV~ 4 / 3 in energy, N 
being the number of particles. This scaling property coincides with that 
for nuclear giant resonances. Contrary to the previous study based on the 
boson-expansion method, the deviation from the harmonic limit is found to 
be almost negligible for Na clusters, the result being consistent with the recent 
experimental observation. 
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Collective phonon excitations are common phenomena in fermionic many-body systems. 
In alkali metal clusters, a very strong dipole plasmon resonances have been observed in the 
electromagnetic response, which are interpreted as collective vibrational excitations of an 
electron gas against a background ions |IJ . These resonances are well described in the random 
phase approximation (RPA) as a particle-hole excitaion, which assumes the harmonic nature 
of the vibrations J| . It may thus be a natural idea to expect the existence of multiple plasmon 
excitations. It is worth mentioning that double phonon excitations have been observed in a 
similar phenomenon in nuclear physics, i.e. giant dipole resonance 

One interesting question is whether the dipole plasmon resonance in metal clusters is 
harmonic enough to allow multiple excitations. Catara et al. used the boson-expansion 
method to discuss two-plasmon excitations in metal clusters ||. They claimed that an- 
harmonic effects are quite large. On the other hand, a comparison between a jellium-RPA 
calculation and the result of time-dependent local density approximation (TDLDA) || sug- 
gests that the anharmonic effects are very small for Na clusters. Recently an experiment was 
performed which addresses the anharmonic properties of plasmon resonances of alkali metal 
clusters [[7]. In this experiment, a large intensity of doubly and triply charged fragments 
were obserbed in the charge distribution of photofragment of Na^ 3 . The ionization energy 
of this cluster lies between the energy of single and double plasmon resonances, and the 
photon energy was set to be slightly larger than the one plasmon energy. The ionization 
was thus energetically possible only if at least two photon are absorbed. The observed ion- 
ization was interpreted as electron emission via multiple plasmon states, suggesting a very 
small anharmonicity, which does not support the prediction by Catara et al. This picture 
was later confirmed theoretically, although a significant correction from a direct two photon 
absorbtion was reported ||. 

The aim of this paper is to show that the anharmonic effects are indeed very small for 
the dipole plasmon resonance in metal clusters, contrary to the prediction by Catara et al. 
To this end, we apply the variational principle for the time-dependent Schrodinger equation. 
The time-dependent variational principle was recently applied to large amplitude collective 
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motions to discuss anharmonic properties of nuclear giant resonances H. Its applicability 



has been tested on a solvable model in Ref. [10]. The time-dependent variational approach 



allows one to estimate relatively easily the energy shift of double phonon state with an 
analytic formula. As we will see, it has the same scaling law concerning the number of 
particle for both nuclear and cluster cases, although the range of the interaction is very 
different for these two systems. 

Consider a system where N electrons interact with each other in a positively charged 
ionic background whose density is given by pi(r). The Hamiltonian for this system is given 
by 

N v 2 1 e 2 N 

where m is the electron mass. V/(r) describes the interaction between the electrons and the 
ionic background. It is given by 

W(r) = _,»/*. "fell, (2) 
J \r — ri\ 

and satisfies the Poisson equation 

VVj(r) = A7re 2 Pl (r). (3) 

As in Ref. [[J, we assume the following wave function to discuss the non-linear dynamics of 
the dipole plasmon resonance of this system. 

Here |^o > is the ground state wave function. The operator Q is the dipole field given by 

N 

Q = Y. z i- (5) 

1=1 

The time evolution of the variables a(t) and (3{t) is determined according to the time- 
dependent variational principle, 

8 [dt < Voflidt - #|* a/3 >= 0. (6) 



This leads to the following two coupled equations || 

$ = -, (7) 
m 

a 2 

ad/3 < VplQlVp > + dp < ^f p \H\^ > +—d < V \[Q, [H,Q))\^ >= 0. (8) 

In order to find the phonon frequency, we requantise these equations of motion with a 
Hamiltonian formulation. One convenient choice of the canonical transformation is |§ 

P^q = P, (9) 
a^p = ma<*p\[Q,[H,Q]}\*p>, (10) 



together with the Hamiltonian of 



H =mr) +U{0) - 

where the inertia and the potential are given by 

M(P)=m 2 <* j3 \[Q,[H,Q]]\* l3 >, (12) 
U(J3) = < * \H\V P >, (13) 

respectively. For the dipole field (|5]), the inertia is easily evaluated to be M(/3) = mN. After 
dropping the constant term, the collective potential is calculated as 



Po being the ground state density. To derive Eq. ([14]), we have used the transformation 

e -mf3lH,Q] {xij ^ z . )e m m ,Q] = ^ ^ z . _ q (15) 

for the dipole field, and the fact that both the kinetic energy and the Coulomb interaction 
among the electrons in the Hamiltonian ([I]) are translational invariant ]nj. Since we are 



interested in the leading order correction to the harmonic limit, we expand the ground state 
density in Eq. fll4]) in terms of (3. Accordingly, we express the collective potential as 
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U(P) = U + ^ 2 + ^P 4 + ---. (16) 

The linear term in the expansion vanishes because of the stability of the ground state, and 
the third order term drops if the spherical symmetry is assumed for the ground state density. 
Using the Poisson equation (§), the coefficients k and k± are evaluated as 

167r 2 r°° 

k = — — e 2 / r 2 drp!(r)p (r), (17) 
6 Jo 

fc4 " 54 6 Jo dr dr ' U8j 

respectively. Here we have assumed that the ionic density pi has the spherical symmetry. 
The expression for k, Eq. (|TTD, was first derived by Brack in the context of the sum rule 
approach 0,0]. 

Equations (ff7|) and (|i~8l) are general expressions and valid for any form of the ionic and 
the electronic density distributions, as long as they are spherical. In order to get a more 
transparent formula, we further simplify them by using the jellium approximation. Here the 
ionic charge density is uniform in a sphere of radius R = r s iV 1//3 , where r s is the Wigner-Seitz 
radius. We thus assume 

Pl (r) = JLe(R-r), (19) 

where 9 is the theta function. Substituting Eq. (|19|) to Eqs. ( |i~7D and ([18]), we obtain 

Ne 2 

k = ^, (20) 



(21) 

r=R 



18 r s dr 

Here we have neglected the effect of spillout of the electron density outside the ionic back- 
ground, which is order of 1/N. To get an analytic expression for k^, we approximate the 
electronic density by the error function as [|13[ 



^^r*(^)- < 22 > 

where a is the surface diffuseness parameter for the electronic density. Substituting this 
density to Eq. we finally obtain 



e 2 N 2/3 



** = -*n=z- ( 23 ) 



3ay/iir 4 s 

We now requantise the collective Hamiltonian Ti. and obtain the phonon frequency. The 
frequency in the harmonic limit is given by 



e 2 



mri 



(24) 



k 

which coincides with the Mie frequency. The leading correction to the harmonic limit is 
given by PJlO] 



"\h, 1 r 2 

E n = ™o - Sf^o 2 = nuj - -^^N^W . (25) 

ofv o\/7T CLC 



Taking the second derivative, the energy shift of the double phonon state is found to be 

= E 2 - 2E X -E = -^Li-N~ A /\ (26) 

Ay/7i ae 2 

Note that the dependence of A&'E/uo on is the same as that for the nuclear giant 
resonances 0. Expressing it as a dependence of the volume of the system iV ~ L 3 , it is 
also the same as that for the photon spectrum in a small cavity |]10[ . All of them scale as 

Let us now estimate numerically the shift of the frequency for Na clusters. Using the 
Wigner-Seits radius of r s = 4 a.u., the Mie frequency (|2~4f) is evaluated as u Q =3.39 eV. 
Bertsch and Eckardt fitted the electronic density obtained by a self-consistent local density 



approximation by Eq. (22]) and obtained a =2.14 a.u. for = 198 | I3| . Assuming that 
the surface diffuseness a has a very weak dependense of N, the energy shift of the double 
phonon state ( p6| ) is estimated as 

A {2) E = -0.45 x N~ 4/3 (el/), (27) 

which is extremely small compared with the phonon frequency Uq in the harmonic limit. 
For example, it is as small as — 1.1 x 10~ 3 eV for A^=92. This result is consistent with the 
recent experimental observation 0] that multiple dipole resonances are easily accessible for 
Na clusters. 
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In summary, we applied the time-dependent variational principle to the dipole plasmon 
resonance of alkali metal clusters to discuss its anharmonic properties. The uniform jellium 
approximation for the ionic density as well as the error function approximation for the 
electronic density lead to a simple analytic expression for the energy shift of the double 
phonon state. We found that the ratio of this quantity to the frequency in the harmonic 
limit scales as iV~ 4 / 3 , which coincides with that for both nuclear giant resonances and for the 
photon spectrum in a small cavity. For Na clusters, we found that the anharmonic effects 
are almost negligible, which is consistent with both the previous TDLDA calculation and 
the recent experimental suggestion. 

The author is grateful to G.F. Bertsch for useful discussions and carefully reading the 
manuscript. This work was supported by the Japan Society for the Promotion of Science 
for Young Scientists. 
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